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ZUR FRAGE DER ADDITIVITÄT DER PARACHORE 
VON LÖSUNGEN 


Von MARIAN PUCHALIK 
Physikalisches Institut der Schlesischen L. Waryúski Medizinischen Akademie, Rokitnica 
(Eingegangen am 10. Januar 1955) 


In der vorliegenden Arbeit wurde die Abhängigkeit der Parachore der Lösungen 
von ihrer Konzentration untersucht. Die Resultate der Messungen zeigen, dass die Addi- 
tivität der Parachore der Lösungen nur in verhältnismässig seltenen Fällen zutrifft. Es 
wurden zwei Gleichungen vorgeschlagen, die die funktionelle Abhângigkeit der Parachore 
der Lósungen von ihrer Konzentration ausdrücken. Bemerkenswert ist die Abnahme 
der Parachore der stark verdünnten Lösungen mit wachsender Konzentration im Fällen, 
in denen das Parachor der gelósten Substanz das des Lósungsmittels weit ubertrifft. 


Das Problem der Additivität der Parachore von Lösungen ist sowohl vom Stand- 
punkt der Theorie wie auch von dem der messtechnischen Praxis wichtig. Die F'eststel- 
lung der funktionellen Abhängigkeit der Parachore der Lösungen von der Konzen- 
tration kann den Ausgangspunkt für die theoretische Deutung der Erscheinungen 
auf der freien Oberfliche der Lósungen bilden. Anderseits ist dieses Problem eng mit 
den Bestrebungen verschiedener Forscher gebunden, das Parachor der gelósten 
Substanz aus den Messungen an Lösungen und an den Lösungsmitteln zu berechnen. 

Das Parachor der Lósungen wird bekanntlich durch folgende Formel definiert 


4 


A 
Pj, = (M, f, + Mofo) Bs (1) 


(Piz bezeichnet das Parachor der Lösung, с» und 013 die Oberflächenspannung, bzw. 
die Dichte der Lösung, M, und M, die Molgewichte, f, und f; die Molbrüche des 
Lósungsmittels und des gelósten Stoffes). 

Zur Berechnung des Parachors der gelósten Substanz bedient man sich der 
Formel von Hammick und Andrew (1929) oder der Formel von Skramowsky (1944. 
nach Soucek 1945). Die erste dieser Formel setzt die Additivitát des Parachors der 
Lósungen voraus: 


Ру» = Р, (1 — fa) + Р», (2) 
(Р, und Р» sind die Parachore des Lósungsmittels und der gelóssten Substanz). 
(379) 
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Nach Skramovsky wird das Parachor der gelösten Substanz aus der folgenden 
Gleichung berechnet: 


(3) 


(Ao ist die Differenz der Oberflächenspannungen der Lösung und der gelösten Substanz, 
До die entsprechende Differenz der Dichten). In der Formel von Skramovsky wird 
es angenomen, dass die Oberflichenspannung der Lósung der folgenden Gleichung 
genügt: 


012 = оу (1 — fe) + о (4) 


Nach den Beobachtungen des Verfassers (Puchalik 1954) treffen diese beiden Voraus- 
setzungen oft nicht zu. Deswegen stellte sich der Verfasser die Aufgabe die funktionelle 
Abhängigkeit der Parachore der Lösungen von der Konzentration dieser Lösungen 
nüher zu untersuchen, mit der Absicht die Grundlagen für die Berechnung der Para- 
chore der gelösten Substanzen auch für solche Fälle zu schaffen, in denen die Addi- 
tivitit nicht mehr zutrifft. 


Das Untersuchungsmaterial 


Zur Lósung des oben angedeuteten Problems bediente sich der Verfasser der 
Resultate seiner vorigen Arbeit (Puchalik 1954) sowie der Messungen von Harkins und 
Grafton (1921). Ausserden wurde speziell für diesen Zweck die Oberflichenspannung 
der Lósungen der a — und f — Naphthole in Aethylaether und in Aethylalkohol 
gemessen. Die genannten Substanzen wurden nach den in früheren Arbeiten des 


Verfasser angewandten Methoden gereinigt (Puchalik 1932, 1954). 


Die Untersuchungsmethode 


Zur Messung der Oberflichenspannung wurde die Methode des. maximalen Luft- 
drucks in Blasen von Cantor-Rehbinder (Rehbinder 1935) angewandt. Als Kapillaren 
wurden bakteriologischen Mikropipetten benutzt, welche sich zu diesem Zwecke als 
gut geeignet erwiesen haben. Als Manometerflüssigkeit wurde raffiniertes Petroleum 
angewandt. Die Änderung des Luftdruckes erfolgte langsam genug, damit man nach 
einiger Übung die Meniskushöhe im Moment des Abreissens-der Luftblase vermittels 
des Komparators leicht ablesen könne. Die Genauigkeit der Messung д betrug с. 
0,3%. Kontrollmessungen gaben gut übereinstimmende Werte mit denen die mittels 
der Steighöhenmethode erhalten wurden. Im Vergleich zur der Steighöhenmethode ist 
die benutzte Rehbindersche Methode genauer, da die Resultate ihrer Anwendung 
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nicht so stark durch die Beschaffenheit der Wände der Kapillaren beeinflusst 


werden. 


Die Konstanthaltung der Temperatur wurde durch ein Thermostat nach Hóppler 
gesichert. 


Besprechung der Resultate 


Die Resultate der Messungen sind in den Tabellen 1 — 8 und in den Dia- 
grammen 1 — 3 zusammengestellt. Aus den gewonnenen Resultaten geht deutlich 
hervor, dass die Additivität der Parachore von Lösungen nur in einigen der untersuch- 
ten Fällen zutrifft. Die untersuchten Lösungen kann man in 3 Gruppen einteilen. 
Zur ersten Gruppe gehóren die Lósungen von Aethylaether in Benzol und von 
Nitrobeuzol in Benzol. Wie aus den Tabellen 1 und 2 ersichtlich ist, kann man die 
funktionelle Abhängigkeit des P,, von f, in diesem Falle durch die vom Verfasser 
vorgeschlagene Formel 


Piz = Р, (1 — fe) + P. (5) 


ausdrücken. Wenn man nämlich aus irgend einer Konzentration der betreffenden 
Lósungen unter Zugrundelegung der von unmittelbaren Messungen stammenden 
Werte für P, den Exponenten Ё berechnet, so bekommt man für die anderen Konzen- 
trationen solche Werte von P,, die mit den unmittelbar gemessenen innerhalb der 
Fehlergrenzen gut übereinstimmen. Damit scheint die Formel (5) für diese Lósungen 
berechtigt zu sein. 


(C242)50-Cg Hg 


Pro = Руди +Pofok 
k *104 


05 1 fe 
Fig. 1 


Anders liegen die Verhältnisse für die zweite Gruppe von Lösungen. Zu dieser 
Gruppe gehören die Lösungen von Aethylalkohol in Wasser (Puchalik 1954) und 
von f — Naphthol in Aethylalkohol. In dieser Gruppe nimmt Р,» mit wachsender 

Konzentrationen anfangs ab, obwohl die Parachore der gelösten Substanzen die des 
Lösungsmittels weit übertreffen. Von einer bestimmten Konzentration ab, steigt Р,» 
linear mit wachsenden f, und bei f, = 1 erreicht den Wert P, für die reine gelöste 
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Substanz (Tab. 7 und Fig. 2). In diesem Falle kann man im Bereiche der nicht allzu 
kleinen Konzentrationen die P,; = 9(f;) — Kurve durch die folgende, ebenfals vom 
Verfasser vorgeschlagene, Gleichung ausdrücken 


Ру = Ko Pi (1 — fa) + Pa V (6) 

wo K, einen echten, konstanten Bruch darstellt. 
Die Parachore der Lósungen von Aethylalkohol in Benzol (Puchalik 1954) und 
von a — Naphthol und f — Naphthol in Aethylaether bilden additive Funktionen der 


C2H50H - H50 


l Ft ME LUE Г 
Fig. 2 


Konzentrationen, Tab. 4 und 5. Diese Additivität der Parachore kann man als einen 
Grenzfall der oben besprochenen Fälle für k=l und k,=1 betrachten. Im Fällen wo 
die Additivität nicht zutrifft ist die Abnahme des Parachors mit wachsender Konzen- 
tration im Bereiche grosser Verdünnungen von besonderem Interesse, obwohl als 
gelóste-Substanz immer die mit dem grösseren Wert des Parachors gewählt wurde. 


42 
53 Ce Hs OH - H0 
52 g= (f. 2) 


Fig. 3 


Ein vortrefliches Beispiel dieses merkwiirdigen Verhaltens der Parachore von 
stark verdünnten Lösungen bilden die wässrigen Lósungen von Phenol. In Tab. 3 
und Fig. 3 sind die Werte der Parachore dieser Lósungen bei verschiedenen Konzen- 
trationen zusamengestellt. Diese Werte wurden aus sehr genauen Bestimmungen der 
Oberflichenspannung von Harkins und Grafton (1921) berechnet. 
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(GH, О — C,H, 


К = 1,042 
fa 


0,0000 


0,1649 
0,2566 
0,3948 
0,6379 
1,0000 


С.Н, МО, — C,H, 


K — 1,093 


0,0000 
0,1821 
0,2545 
0,3935 
0,6074 
0,7603 
1,0000 


Tab. 1. 
t = 20°C+ 0,1°C 
Py Р, 
— —ͤ—. S 
206,0 
204,1 210,8 
204,3 211 
205,3 212 
206,0 210 
210,8 210,8 
Таб, 2- 
t = 20°С + 0,1°С 
206,0 
210,3 262,0 
217,9 269 
222,5 263 
232,6 254 
246,6 260 
262,0 262,0 
Tab: 3: 


C,H,OH — H,O nach Harkins und Grafton 
t = 20°С 


0,00018 


0,00090 
0,00181 
0,00366 
0,00541 
0,00715 
0,00892 


0,9987 
1,0002 „ 
1,0009 ,, 


Tab. 4. 
a — C,,H,OH — (C,H,),0 
Р, = 209,3 


t = 20°C + 0,1°C 


Tab. 5. 
В — Co H, OH — (C,H,),0 t= 20°С + 0,1°C 
PI = 200, 3 


fa 


0,0602 


0,0936 
0,1070 
0,1277 


Tab. 6. 
a — Cy) H,OH — C,H,0H 
P, = 126,9 


t = 20°С + 0,1°С 


0,0056 


0,0112 127,4 
0,0230 129,5 
0,0601 138,0 


Tab. 7. 


В — CoH,0H — C. HOH : = 20°C + 0,1°C 


Dynen 


0,0052 20,04 

сш 
0,0180 теба 5 123,9 
0,0262 1819 125,5 
0,0394 18,03 „ 127,6 


Tab. 8. 
t = 20°С + 0,1°C 


а— СьН;ОН — C,H, 
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Leider kann man wegen der schweren Löslichkeit des Phenols in Wasser in die- 
sem Falle den Gang der Piz = 9(f,) — Kurve auch für grössere Konzentrationen 
nicht bestimmen. 

Dasselbe betrifft die Lósungen der a — und f — Naphthole in Wasser und in 
Benzol. 

Da aber die Werte der Parachore verschiedener Lósungen im Bereiche sehr 
kleiner Konzentrationen im Allgemeinen nicht bestimmt sind und sie genügen auch 
nicht den Formeln für gróssere Konzentrationen, wird vom Verfasser das Problem der 
genaueren Untersuchung solcher Lösungen wie auch das der theoretischen Deutung 
dieser Erscheinungen gestellt. Bezugnehmend auf die in Tab. 6 dargestellten Ergebnisse, 
muss die Sonderstellung der Parachorwerte der Lósung von a — Napthol in Aethylal- 
kohol näher besprochen werden. In diesem Falle bildet das Parachor der Lösung keine 
streng additive Funktion der Konzentrationen, es nimmt aber monoton mit wachsen- 
der Konzentration zu. Bemerkenswert ist, dass die Oberflächenspannung in diesem 
Falle innerhalb der Fehlergrenzen als konstant anzusehen ist. 

Die zur Untersuchung nótigen a — und 8 — Naphthole wurden dem Verfasser 
durch das Institut für Physiologische Chemie der Schlesischen Medizinischen Aka- 
demie geliefert. Dem Leiter dieses Instituts, Herrn Prof. Dr Stanislaus Józkiewicz, 
spricht der Verfasser seinen verbindlichsten Dank aus. 


КРАТКОЕ СОДЕРЖАНИЕ 
Пухалик, К ВОПРОСУ OB АДДИВИТЕТЕ ПАРАХОРОВ РАСТВОРОВ. 


В настоящей работе исследована зависимость парахоров растворов от их кон- 
центрации. Результаты измерений указывают на то, что аддивитивность пара- 
хоров растворов сравнительно редко случается. Были предложены два уравнения, 
которые выражают функциональную зависимость парахоров растворов от их 
концентрации. Замечательное уменьшение парахоров сильно разбавленных рас- 
творов с возростающей концентрацией в случаях, в которых парохор разбалвен- 
ного вещества много превышает парахор растворителя. 
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AN INTEGRATING 
APPARATUS FOR THE SCHRÓDINGER EQUATION. II 


Ву К. Antonowicz 
Department of Physics, Nicholas Copernicus University, Тогай 


(Received February 17, 1955) 


Several improvements have been introduced in the integrating apparatus for the 
Schrödinger equation, (Antonowicz 1953). The apparatus was calibrated by means of the 
Morse potential and its harmonic approximation. For verification purposes eigenvalues 
were determined for the harmonic oscillator and the Morse potential. The eigenvalues 
obtained by means of the apparatus agree within 1 per cent with the theoretical values. 
One of the eigenfunctions for Morse's potential was compared with a theoretical curve. 
А very good agreement between them was found. 


I. Introduction. In an earlier paper (Antonowicz 1953) a description of an inte- 
grating apparatus for the Schródinger equation was given. The apparatus is based on 
the analogy between the Schródinger equation and that of the motion of a magnetic 
needle in a homogeneous magnetic field, i.e. on the analogy between the equations 


= “е Dt (1) 
r H, Hija-o @ 


After our first paper appeared, we had the possibility of becoming acquainted with 
a paper of Bullard and Moon (1931), from which we learned that the principle on which 
the construction of our apparatus was based was similar that already used by these 
authors. However, the realization of the principle was different. Apart from technicali- 
ties (e.g. the use of a movable needle instead of a movable coil) there were differences 
resulting from the desire of eliminating as much as possible the share of the experimen- 
ter in the process of solving the equation. 

By means of the described apparatus we were able to obtain eigenfunctions 
(of the harmonic oscillator) with a very great accuracy. We met, however, with serious 
difficulties in determining the eigenvalues and calibrating the apparatus. The stability 
was too bad and the sensitivity for external disturbances too great to enable accurate 


(385) 
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measurements to be made. Thus necessary steps were taken to improve the stability 
and precision of the apparatus; namely 

(1) by application of an astatic needle system; 

(2) by stabilizing the D. С. amplifier and the speed of motor rotations; 

(3) by construction of an apparatus for obtaining accurate ,potential curves". 

A short description of the technical improvements, the method of calibrating 
the apparatus, as well as the results obtained by integrating the Schródinger equation 
for Morse's potential are given below. 

II. The magnetic needle. The oscillations of a magnetic needle suspended ona 
cocon fibre are strongly damped if the amplitudes are rather great, whereas at small 
amplitudes external disturbances cause continual oscillations of the needle. The 
external disturbances are of two kinds: mechanical and magnetic. The mechanical 
disturbances might be eliminated by a shock-proof suspension system. It is better, 
however, to eliminate them by using a moving system with such a high moment of 
inertia that the system becomes insensitive to these external mechanical disturbances; 
an additional advantage of such a system is that the damping of oscillations becomes 
comparatively small and the analogy between equations (1) and (2) closer. In our 
apparatus the moment of inertia of the astatic needle system has been increased by 
means of a lead cylinder (external diameter 16 mm, inner diameter 6 mm and 5 mm 
in height). 

Magnetic disturbances are eliminated by using an astatic needle system, the needles 
being suspended one above the other at a distance of 10 cm. The needles of 2x 2x8 mm 
were made of magnetic steel. Each needle is in the homogeneous part of the magnetic 
field given by two Helmholtz systems of coils. Each of the Helmholtz systems has 3 
independent windings. The róle of the first winding is to give H(t), the second and the 
third ones are used for precise adjustment of H, for finding the eigenvalues. We are unable 
to measure currents by means of our milliammeters with an exceeding precision 1 per 
cent.; thus it is very difficult to adjust the currents with a better precision. In searching 
for the eigenvalues, however, it is necessary to regulate the intensity of the field with 
a precision of at least 0,01 per cent. This difficulty is solved by the use of two windings. 
One of them with a great number of turns is situated near the needle, the other one 
having but a few turns is placed at the farthest end of the coils. The same current in 
the first winding produces a considerably stronger field than it would produce flowing 
in the second one. Searching for the eigenvalues, we first regulate the intensity of the 
current in the first winding, with as much precision as the milliammeter permits, 
and then in the second one, which produces small changes in the magnetic field even 
if changes of the current are considerable. In this way a subtle regulation of the magnet- 
ic field is achieved. The measuring of the magnetic field is carried out independently 
of the measuring of the current (see below). 

III. The time-dependent magnetic field. The field H(t) is produced by the ampli- 
fied current of a photocell. The intensity of light falling on the photocell is regulated 
by a disc (, potential curve“) revolving at constant speed, the edge of which is cut in 
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the form of the U(x) function in the polar system. To make sure that all the discs are 
made on the same scale and with the greatest precision possible, a special apparatus 
has been constructed, (see Fig. 2). 

It consists of an iron plate, 50 cm in diameter, that may be turned round axis 
O with an angular scale at the edge divided into 0,5?. Above the plate a micrometer 


Fig. 1. Astatic needle system. DCA — direct current amplifier. A — resistor for coarse regulating of Ho. 
C — resistor for fine regulation of H,. B — switch used to stop the needle. D — switch used to switch 
off H (t) 


screw is fastened which permits to shift slide B by means of knob C along the radius 
of the plate; the position of the slider can be read with a precision not worse than 0,01 
mm. The slider is provided with a puncher. When puncher B is pressed it makes a small 
trace (hole) on the paper placed on the plate. The equation of the curve being known, 
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the calculated points may be very accurately put down on paper and thus enable a very 
accurate „potential curve“ to be cut out. 

The direct current amplifier must be characterized by an exact reproduction of 
currents and, apart from this, by a great stability and considerable current output. 
In order to ensure stability of the amplifier its power sources have been stabilized. The 
incandescent tubes were fed by a battery of accumulators. The plate voltage were fed 


Fig. 2. The apparatus for marking the potential curve. O — axis of the plate. А — micrometer screw. 
B — puncher. € — knob of the screw 


by a rectifier stabilized by a tubes system which, in turn, was fed by mains through 
a magnetic stabilizer. Our amplifying system is generally known as а „bridge circuit“. 
A certain shifting of the working point is observed which, however, does not affect 
the precision of measurements, since a suitable correction in the eigenvalue can be 
easily made. 

The stability of the speed of the motor plays an important rôle. A small converter 
fed by direct current from a battery of accumulators is used as motor. The dynamo of 
the converter feeds a stabilivolt. As it is known, small voltage changes produce great 
current changes in the stabilivolt. Thus the load of the converter is changing and regu- 
lates the motor speed. The speed is checked by means of a milliammeter introduced 
into the converter circuit. 4 0,01 per cent. change in speed causes c. 1 mA change of 
the current. А resistance on the side of motor feeding permits to regulate the motor 
speed up to the required value. ^ 

IV. The eigenvalues. The value of the constant magnetic field H, by which the 
function satisfies the boundary conditions, determines the eigenvalue. H, is not directly 
proportional to the eigenvalue; it must be corrected by subtracting the field value 
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due to the constant component of the amplifier current, the turning moment of the 
cocon fibre, as well as by the value of the earth magnetic field with which the movable 
system interacts in spite of the astaticity of the needle system. Thus a certain value Hy, 
must be subtracted from Ho in order to obtain a value proportional to the eigenvalue E,. 
The value Ho is easy to determine in the case of the harmonic oscillator, in which 
case the difference between the energy levels is constant and equal to Av, the lowest 
level having the value /v/2. Thus subtracting from the lowest measured eigenvalue 
the value corresponding to Av[2 yields Ни. The value Hy, is then subtracted also 
from all the other eigenvalues measured. 

The magnetic field corresponding to the eigenvalue is determined in the following 
way: we let the needle oscillate in the field Ну; in this case the curve obtained by 
means of the apparatus is a sine curve, from which we calculate H,. Table I contains 
the eigenvalues thus determined for the harmonic oscillator in the units hy = 1. 


Table 1. Eigenvalues of the harmonic oscillator obtained by means of the apparatus 


0,5 0,500 


Quantum number 


0,5 per cent 
8 
3 207 
Об 
Ole: a. 
02x 5 
Ole, 
Ol 


© -з с C14» м س‎ о> 


V. Calibration of the apparatus. In order to calibrate the apparatus eigenfunc - 
tions and eigenvalues have been found for the Morse potential (Morse 1929) as well 
as for its harmonic approximation. Аз it is known this potential has the form 


U(r) = D(1 — e)2 


f 


and the frequancy of small classical oscilations about the point r, is 


_ а (2D\" 
MET a 


Therefore, if we take then constant factor in the expression for the potential energy 
of the harmonic oscillator as equal to 
2л? veli = aD, 


such a harmonic oscillator will represent the zero approximation to the Morse poten- 
tial. ў 
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We assume the following values of the parameters: 


D = 50 000 cm 
а ШЕ IN AU. 
e eet n 


The potential curves were made by means of the special apparatus described above 
(Fig. 2) in the scale 1º = 0,01A. U., and 1 mm = 1,10? стг! = 9,9275 1079 erg. 
The eigenvalues obtained for the potential curves were in arbitrary units. Since the 
difference between the eigenvalues of the harmonic oscillator is equal to Л», and the 
distribution of the eigenvalues for Morse's potential is given by 


2% 2 2 
By = — D+ nei) Baie Yo (n+ 2) 


4D 


E, and h» being given in arbitrary units, it was possible to determine immediately 
121214 in the same units. The absolute values of Ay and h?9¢/4D may be also de- 
termined in absolute units for a reduced mass u, which results from the ratio of the 
magnetic moment of the needle m to its moment of inertia B, i. e. from m/B (Eq. (2)). 
Thus it is possible to solve the problem for a certain mass и which may be determined 
from the above data. We have obtained by means of the apparatus: 


hv, = 1129 in arbitrary units 
һ2у2/4) = 20,7 in the same units. The ratio of these values does пої depend on the 
units and is equal to 54, 541. 

Exact solutions of the Schródinger equation for the Morse potential are. 


a?D\* 
im ng) 


4D 8r u 
Hence 
hw V 322 uD 
ar s 
4D 


and after substituting the assumed parameter values and the values of universal con- 
stants, we obtain u = 1,20 10-24 .g. 

It is obvious that, the mass corresponding to the ratio m/B being known, it is 
easy to adapt the apparatus to another mass by changing the value m/B (practically 
by changing the moment of inertia B). The same effect may be reached without chang- 
ing this moment of inertia, either by changing the scale of the potential curve, or by 
changing the speed of disc rotations. Thus the apparatus can be adapted to the solution 
of problems with any mass. 

Assuming the mass u = 1,2 10-24g, we can calculate all the quantities we are 
interested in and compare them with those given by the apparatus, namely 
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hyo — in arbitrary units, obtained by means of the apparatus = 1129; 

hy, — calculated theoretically = 7,285 10-13 erg, theorefore the arbitrary energy 
unit = 6,4526 10-18 erg. We can now calculate the eigenvalues for the assumed Morse 
potential, e.g. in ergs, both from the data obtained by means of the apparatus and 
from Morse's expression for eigenvalues. The results are given in Table II. 


Table. II. Eigenvalues for Morses potential 


deeem Eigenvalues obtained by Eigenvalues calculated Di oħtna 
number we “р ГК فی‎ and calculated eigenvalues 

0 —95,728 — 95,666 0,08 рег сепї 

1 —88,400 —88,649 QS у gy 

2 — 81,950 —81,949 D 

3 — 75.655 — 75,417 QS > 

4 —69,405 —69,202 0,4 „ „ 

5 — 63,324 — 563,255 о = 

6 —57,906 —57,576 MW в = 

U —51,751 — 52,164 ОВИЕ 


The absolute error is nearly the same for all eigenvalues, therefore the relative error 
for small eigenvalues is rather greater; the relative error in the differences between 
neighbouring energy levels reaches several per cent. 

VI. Eigenfunctions for Morse's potential. In the previous paper (1953) eigen- 
functions of the harmonic oscillator, obtained by means of the apparatus, were publ- 


Fig. 3. Eigenfunction obtained by means of the apparatus. The points inside the circles mean theoret- 
ical points. The dotted line is the lower part of the (Morse's) potential curve 
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ished and showed very good agreement with the points calculated theoretically. However, 
the shape of the eigenfunctions of the harmonic oscillator does not depend on the 
specific values of the constants in the potential functions. It was therefore interesi 
ing to verify the accuracy of the eigenfunctions of Morse’s potential, whose shape is 


Nod T OEE 


Fig. 4. 8 eigenfunctions of Morse’s potential. All the curves are redrawn for reproduction purposes 


determined by the constants. The eigenfunctions given by the apparatus are not 
normalized. Thus we are free to choose arbitrary one constant, e. g. the normal- . 
izationfactor. The radial eigenfunction for Morse’s potential is given by the 
expression: 


-de а (r- ro) en K—2n-1 -a(r—ro) 
у= С.е lde f 8 © рк—зя-1 [о e-, 


п—1 
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where in the case in question 


qo 20N2uD _ n 
ah 

а= 17 

Ty = 1,3 

* = 2d = 54,2 


For verification purposes an eigenfunction of three zeros was chosen, i. e. for n = 3. 
The dimensions of the curve obtained directly from the apparatus are too small to 
make accurate determination of its errors possible. Therefore the curve was magnified. 
The length of the magnified curve was 1 m and the maximum amplitude 0,15 m. 

By substituting the values of the constants, we obtained the following expression 
for the eigenfunction 

y = 1.1015 e7271 61 6-53) . [5-1 671,323 . [15.922 9-51 671.9 
— 44,241 e. 719 4 40,159 e- L7 0759 — 11,9049] mm 
It is convenient to put 
„„ 
to get a simpler expression, namely 
y = 1.1015 e271 4238 [15,922 х3 — 44,241 x? + 40,1596 x — 11,9049] 

In spite of this simplification the calculation of the points of the curve is very 
tedious. Therefore the calculations were carried out only for == 0,9; 1,0; .. I, 9 
A. U. The teoretical points appeared to be situated within the breadth of the line 
(Fig. 3). The points inside the circles means theoretical points. 

Fig. 4 shows Ө eigenfunctions of Morse's potential, from which the fourth was used for 
checking the exactitude of the apparatus. 

The author expresses his gratidude to professor A. Jablonski for helpful discus- 
sions and support throughout the work. 


KPATKOE СОДЕРЖАНИЕ Е 
Антонович. ПРИБОР ДЛЯ ИНТЕГРИРОВАНИЯ УРАВНЕНИЯ ШРЕ- 
ДИНГЕРА II. 
Введено много улучшений в приборе для интегрирования уравнения Шре- 
дингера. Прибор был калиброван при помощи потенциала Морса и его гармо- 


` нического приближения. Для проверки действия прибора были определены зна- 


чения собственные для гармонического осциллятора и потенциала Морса. Значе- 
чения собственные полученные при помощи прибора, согласуются с теоретиче- 
скими значениями с точностью до 10%. Одна из функций для потенциала Морса 
была сопоставлена в теоретической кривой. Получена очень хорошая согласо- 
ванность между обеими кривыми. — 
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DEPENDENCE OF THE DIELECTRIC PROPERTIES OF CERAMIC 
BaTiO, FOR HIGH FREQUENCY CURRENTS ON THE TECHNOLO- 
GY OF THE PREPARATION OF SAMPLES 


By M. Jeżewski AND T. PiecH 
Institute of Technical Physics, Mining Academy, Cracow 
(Received April 4, 1955) 


The dependence of the dielectric constant and the ceramic BaTiO, losses on the 
way of preparing samples has been investigated. The method of measurements of the 
dielectric constant of materials having very high dielectric constants has been worked 
out with the help of high frequency currents, as well as the method of measurements of 
losses. Four groups of BaTiO, samples, marked 4, В, С, D, have been prepared, each 
of them being obtained from the previous one by means of thorough grinding, pressing, 
and resintering at a temperature of about 1350?C. The last group D remained particularly 
long in the oven: about 130 hours, including the time of heating and cooling. It has been 
proved that the dielectric constant of each successive group of samples is at all tempe- 
ratures considerably higher than the dielectric constant of the previous group. The losses 
of samples in the row A, В, C, D decrease gradually but not very significantly. 


Introduction 


The examination of the dielectric properties of ferroelectrics of the BaTi0, type 
are usually carried out with ceramic-materials, which can be formed according to the 
desired shape and size, which is particularly important in practical applications. The 
comparison of the results of investigations carried out with single crystals and ceramic 
bodies shows that materials having multicrystalline structure behave, in general, 
similarly to crystals (Merz 1949, Rzhanov 1949, Mason and Matthias, 1954). For 
instance the dependence of the dielectric constant of the multicrystalline BaTiO; on 
different parameters is the same as for single crystals, and the difference exists only 
in the absolute value of this constant. This is due to the fact that while examining the 
ceramic samples one measures only an average value of the constant of the crystal- 
lites, their axis directions being distributed at random without any privileged direc- 
tion in space. Yet in reviewing the already comprehensive literature on this subject, 
it can be stated that the data of various authors concerning the dielectric properties 
of the ceramic BaTiO, in weak fields, rather considerably differ in values (Roberts 
1949, Rzhanov 1949, von Hippel 1950, Pajak 1950, Kazarnovskij 1954), in spite 
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of the fact that the samples were prepared by similar methods and that the dielectric 
constant of this material has been found to be independent of the current frequency 
up to f = 108 Hz. (Powles 1948, von Hippel 1950, Piekara 1950). It seems that the 
cause of these divergencies can by no means be ascribed to the possibility of the 
samples, examined by various authors, being prepared from different initial materials, 
thus containing possibly some significant impurities. Although it is well known that 
materials having a high dielectric constant are sensitive to the presence of certain 
additions, yet there is little probability that the commercially pure components, 
used for the synthesis, might have contained impurities in amounts able to cause 
such considerable divergences. Investigations were also made (Graf 1951) proving 
that there exist only small differences between the properties of the samples obtained 
from the initial materials with a high degree of spectroscopically examined purity 
and those of the samples prepared from commercially pure materials. 

It is, however, possible that the dielectric properties of the ceramic BaTiO, de- 
pend on the technique of preparation of the samples. The aim of the present paper 
is to ascertain wheter these properties depend on the way of preparing and the time of 
sintering the samples, and if they do, in what a manner. In order to find it out the 
dielectric constants and the dielectric losses of the ceramic samples prepared in va- 
rious ways (described below) from the same initial materials have been measured at 
different temperatures. Measurements were made by means of a 5 MHz frequency 
current at temperatures ranging from room temperature up to temperatures conside- 
rably higher than the Curie point (of about 170?C). 


Technology of preparing the ceramic samples 


Barium carbonate (Ba'TiO,) indicated as „chemically pure“ as well as titanium 
oxide (TiO,) were weighed in the molar ratio of 1 : 1 and mixed thoroughly by long 
grinding in a porcelain mortar. From this mixture tablets 14 mm in diameter and 
about 5 mm thick were formed under a pressure of 2400 kG/cm?. The sintering of 
these tablets took place in an electric furnace of fire clay at a temperature of 1350?C. 
Heating up to this temperature lasted about 4 hours and for the same period of time 
the samples remained in the oven at peak temperature, while cooling lasted about 
10 hours. According to data from literature (Trzebiatowski 1952, Wojciechowska and 
Damm 1952) in these conditions the synthesis of barium titanate should be 
completely accomplished. In order to avoid diffusion of foreign atoms from the liner, 
on which the samples were placed in the oven during the sintering, a liner from the 
same material in the shape of a brick 6 x 6 x 1 cm in size has been first prepared. 
This liner was sintered in the same way as the tablets on a nickel liner, heated befo- 
rehand in order to bring about the oxidation of the surface. All the samples were 
sintered in the oven on this liner. Ready ceramic samples were then ground on a buff 
wheel to the required thickness. 
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Samples sintered once im the above described manner will be further called 
samples 4. Besides, samples sintered twice (B), three times (С), and four times (D) 
were prepared in the following manner. Samples 4 have been pulverized in a porce- 
lain mortar and the powder thus obtained has been sifted through a sieve of 10000 
eyes/cm?. From this powder tablets were formed under the same pressure as before 
and resintered under the same conditions. In a similar way were prepared 
the samples C out of the B ones, and the samples D out of the tablets C, with 
the difference that the sintering of the latter took place in slightly changed condi- 
tions. Namely, samples D remained in the oven for 130 hours, while the period of 
heating up to the highest temperature (1340°С) amounted to 72 hours. The samples 
remained 8 hours at peak temperature and the cooling process lasted 50 hours. The 
sintering was performed in an electric furnace of fire clay in the Fireproof Material 
Factory in Skawina. 

The samples 4 and B differed distinctly at first sight. Samples 4 were light and 
brittle and microscopic observation showed clearly their porous structure. The sam- 
ples B having the same size were much heavier and after being ground had a smooth 
surface. They had also a greater hardness. The more times the samples were sintered 
the harder they were and the more compact structure they possessed. Ceramic sam- 
ples are, however, always a more or less porous material. Аз far as its dielectric pro- 
perties are concerned the porous material can be regarded as a system of small par-, 
ticles of the given material, immersed in a medium of a different dielectric constant — 
the air. The dielectric constant of such a system will depend on many factors, such 
as the dielectric constant of the components, the shape of the „pores“ filled with 
air, and the ratio of the volume of the air contained in the pores to the volume of 
the whole sample. Only the last of these factors can be established by direct measure- 
ments. Its measure is the so called „absolute porosity“ determined by means of the 
following formula (Polish Ceramic Norms, 1936) ; 


P% = | u Э 100 


where d, = m/v is equal to the ratio of the mass of the dry sample (т) to the volume 
of the sample (v), determined by means of a hydrostatic balance after its previous 
saturation with water, and d is the density of the material out of which the sample was 
made, determined picnometrically after its pulverization. 

Porosity measurements of the examined samples gave the following results: 


Samples Р% 
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Method of measurements of the dielectric constant and of the dielectric losses 


The measurement of capacity was made by means of a specially elaborated reso- 
nance method, which permitted to carry out the capacity measurements with a great 
accuracy even at high losses of the measured condensers, as well as to determine the 
dielectric losses. 

A diagram of the equipment is shown in Fig. 1. The crystal oscillator I of fre- 
quency 5 MHz acts by induction on circuit II set up for resonance. The coupling bet- 


to mut.amm 


Fig. l. Apparatus for measuring the dielectric constant. I — crystal oscillator, producing vibrations of 
frequency 5 MHz, II — resonance circuit, III — resonance amplifier, IV — rectifying circuit connected 
to the milliamperemeter 


ween these circuits was extremely weak in order that the resonance curve be as little 
deformed as possible, and that the voltages in the resonance circuit be as low as 
possible. It was proved that at tuning the resonance circuit does not influence signi- 
ficantly the intensity of oscillations in the oscillator circuit. The coupling between 
the circuits was purely magnetic. All parts of the circuit which could have appreciable 
variable voltages have been taken away as far as possible and carefully screened. The 
voltage on the coatings of the variable condenser in the resonance circuit measured 
with the help of a valve voltmeter reached, at resonance, the value of 4,5 V, when 
there were no losses in the circuit. 

As oscillations in the resonance circuit must have been very weak, a single lamp 
resonance amplifier (III) was used for determining the resonance. The resistance 
circuit acted purely by induction on the grid cathode circuit of the resonance ampli- 
fier, whereas the circuit which is in the anode circuit of the amplifier, tuned to the 
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generator frequency, acted purely by induction on the aperiodic circuit IV, which 
contained a rectifier lamp and a milliamperemeter. Individual circuits interacted on 
one another only by means of coupled coils. In order to avoid the action of individual 
circuits upon other circuits, the whole equipment has been mounted in an iron-plate 
box with double walls, divided also by means of iron- plate walls into separate cham- 
bers. The walls have been marked on Fig. 1 with a dash-line. Outside the box there 
was only a knob rotating the variable condenser in the resonance circuit, a milliampe- 
remeter and the lamps. A disc with a scale has been placed on the axis of the variable 
condenser. Through a small window in the upper wall of the box the rotation angles 
could be read. А nonius placed close to the scale enabled the reading of the tenths 
parts of a scale division. Particular investigations showed that each of the circuits 
acted only on that circuit with which it was directly coupled. 

The measurement consisted (1) in reading the two capacities of the variable con- 
denser with disconnected test condenser, while the deflexions at both sides of the 
resonance point are equal, (2) the same readings after the test condenser has been 
switched in the circuit, while the deflexions of the amperemeter are the same as 
before. Subtracting from the sum of the two former figures the sum of the two latter, 
we obtain the double value of the condenser capacity in degrees of the scale of the 
variable condenser. The result obtained is independent of the conductivity and of the 
losses of the test condenser (Jezewski 1927, 1933). The variable condenser possessed 
a capacity varying from the initial capacity C, to С, + 100 pF. As the capacity of 
the test condenser with BaTiO, samples exceeded, near the Curie point, the value of 
100 pF, the scale of the variable condenser has been raised by connecting to it in paral- 
lel two constant condensers: 50 pF and 100 pF. Owing to these additional condensers 
it was possible to measure capacities up to 250 pF. The additional condensers were 
switched on in the circuit by means of pin-switches, excluding any capacity variations 
due to displacement of any parts with respect to one another. 

When a condenser with losses has been switched on, the resonance curve was 
lowered. Losses were established by measuring the resistance connected in parallel 
to the variable condenser, which caused identical lowering of the maximum of the 
resonance curve. On account of the small sensitiveness of the resonance curve to 
resistances of about 1 megohm, this method, inaccurate at low losses, becomes more 
accurate at higher ones. In the case of resistances of, say, 1 megohm, the accuracy of 
measurements amounted to about 20%, at resistances of, say, 10 КО (which 
corresponds with the losses of the examined samples near the Curie point) it 
reached 1%. | 

As the absolute value of the dielectric constant had to be exactly measured, 
and not only its relative changes, a plate condenser with three electrodes with its 
outer coatings earthed has been used for the measurements. In such a condenser the 
inner coating is well screened by the outer ones, due to which the influence of the 
stray capacities in respect to the surrounding space is insignificant (Jellonek 1953). 
For calculating the capacity of the condenser filled with air, the following formula 
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has been employed (Kohlrausch, Praktische Physik), 


72 r d 


where r is the radius of the circular inner coating, а — the thickness of the coatings, 
d — the thickness of the inner electrode (k). The capacity measurements could be 
carried out with an accuracy of 0,1 pF. Yet on account of the small accuracy in the 
calculation of the capacity of the condenser filled with air, which is dependent on the 
diameter of a very small inner coating as well as on the 
thickness of the tablets, one can hardly ascribe to the 
values of the dielectric constant a greater accuracy 
than 195. 

The technical details of the condenser used, are 
shown in Fig. 2. Metal plates 4 and B pressed together 
two BaTiO, tablets (C), of equal size, silver-plated on 
their whole surfaces adjacent to the plates, while on their 
inner side there were silver electrodes having the shape of 
circles of about 0,6 mm in diameter. Between these elec- 
trodes there was a brass disc having the same diameter 
as the electrodes. Plates 4 and B were connected to the 

с earthed coating of the variable condenser. The electrode 
Fig. 2. Condenser for measu- : 4 ams 
ring the dieset constant of k could be inserted in the circuit by means of a metal 
the ceramic BaTiO, Plates А pin with a very thin edge. This rod could be raised and 
and B pressing the sample C, kept at constant distance from the inner coating by an 
connected to the earthed coa- adequate stopper (z). The tablets were silvered by means 


Ung of the variable condenser e firing samples covered with a silvering paste in sui- 
in the resonance circuit, k — 
table places. 


inner electrode, z — stopper 
of the pin, switching on the Making tiny inner electrodes presented a particular 
test condenser in the circuit difficulty. On account of the necessity of dealing with 
small capacities and simultaneously with regard to a very 
large value of the dielectric constant of the examined material, it was necessary to re- 
duce the diameter of the inner coating below 1 mm. It was very difficult to silver 
such a small coating, which should have had exactly the shape of a regular circle. 
These technical difficulties have been surmounted in the following manner. By means. 
of a cutter of a suitable diameter a cylindrical cavity about 0,2 mm deep was drilled 
in a prepared sample, the inside of the cavity was carefully silvered, fired in an oven 
and then the surface of the plate was slightly ground in order to remove the traces of 
silvering around the cavity. The coating was left a little sunk (0,1 mm); which enabled 
later on the central adjustment of the electrode k. The diameter of the electrodes was 
measured under the miscroscope by means of an eyepiece micrometer. 
The test condenser was placed in an adequately built electric oven, whose tempe- 
rature could be changed from room temperature up to temperatures above 200°C. 
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The temperature of the condenser was measured by means of a gradauted thermo- 
couple, one junction of which touched the surface of the tablet. The second junction 
was held in ice. The temperature was measured with an accuracy up to 0,1°. 

At such a high current frequency a marked influence of the inductance of the 
connecting wires upon the capacity measurements could be expected, the results of 
the measurements carried out by the resonance method were, therefore, compared 
with the results obtained by the accurate bridge method at a frequency of 4.5 kHz. 
The results of these measurements proved to be compatible. Thus it has been esta- 
blished that the influence of the inductance of the wires and condensers is smaller 
than the experimental errors. Evidently the inductances of the variable condenser 
in the resonance circuit together with the connecting wires, and the inductances of 
the test condenser compensated one another. 


Results of the measurements 


The results of the measurements of the dielectric constant as a function of the 
temperature were scheduled in Table I and in Fig. 3, and the results of the dielectric 
losses in Table II and in Fig. 4. As the differences in values of the losses for samples 
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Fig. 3. Dieletric constants of samples 4, B, and D of ceramic BaTiO, as functions of the temperature. 
Crosses indicate the values of the dielectric constant of the ceramic BaTiO, according to the mea- 
É. surements of Roberts 
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A, B and C lie within the limits of error, the values of the losses for samples В and С 
were not given here. Also for the sake of transparency of the drawing, the Pau 
representing the values of the dielectric constant of sample С, which do not differ 
considerably from the corresponding values for samples B, were not shown in Fig. 3. 


20р 509 100° 150° 


Fig. 4. tan д for samples 4 and D of the ceramic BaTiO, as functions of the temperature. Crosses indicate 
the losses according to the measurements of Roberts 


Table I 

A -B D 

tC °С tec 
21,0 129 21,0 1190 - 20,0 1475 
44,5 120 51,5 1104 58,0 1280 
50,0 119 60,0 1108 68,7 1285 
57,1 119 71,25 1124 85,0 1446 
69,5 120 89,0 1184 101,5 1661 
81,0 123 100,6 1353 110,0 1895 
97,5 138 111,0 1850 115,0 3615 
103,5 146 1147 3014 116,0 5780 
109,0 147 116,0 3643 117,0 6920 
112,0 188 117,0 3686 117,5 7051 
116,0 318 119,0 3923 126,0 6100 
119,2 372 120,0 3918 128,7 5975 
120,5 365 121,5 3860 139,7 4240 
122,0 349 124,5 3725 152,5 3015 
- 1340 272 125,0 3680 154,0 2980 
146,7 218 145,0 2502 162,2 2430 


174,5 154 173,5 1612 172,2 2045 


"a. 
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Table II. 
A D 

t°C tan. Ó t°C tan. д 

18,0 0,052 22,0 0,049 

33,0 0,047 36,0 0,035 

55,5 0,036 45,0 0,038 

76,0 0,043 72,0 0,028 

87,5 0,059 103,5 0,069 
107,0 0,094 115:5 0,087 
118,0 0,092 116,0 0,085 
127,5 0,079 125,5 0,073 
137,5 0,068 153,0 0,044 

Conclusions 


From the figures above, as well as from the Figs. 3 and 4 it is easy to prove that 
the way of preparing and the time of sintering the samples has, in general, no 
influence on the character of the dependence on temperature of the dielectric 
constant and the losses of the ceramic BaTiO,. According to the results of investi- 
gations of the majority of authors, with increasing temperature, the dielectric 
constant, decreases at first from room temperature upwards, reaching a flat minimum 
at a temperature of about 60°С, then, from about 110°С it increases rapidly, reaching 
a maximum at the Curie point, which for all examined samples possesses the same 
value (within the limits of accuracy of measurements in the interval of 119,5°—120°C). 
The values of the dielectric constants of all samples decrease above the Curie point, 
but the slope of the curve is much smaller than below the Curie point. For compa- 
rison the results of the measurements of Roberts (1949) have been marked with cros- 
ses in Fig. 3, and it can be seen that they usually fall between the results for our sam- 
ples B and D. | j 

As can be seen from Fig. 4, the losses show a somewhat similar temperature 
dependence. Starting from room temperature, the losses decrease first with increasing 


temperature, reaching a minimum between 60° and 70°C, and later on they increase, 


reaching a maximum below the Curie point. For various samples this maximum lies 
at somewhat different temperatures. From Fig. 4 it can be seen, however, that this 
maximum comes nearer the Curie point for samples having a more compact structure 
(sample D). The influence of the way of preparing and the time of sintering the sam- 
ples on the values of losses is rather small. Samples sintered many times do not show 
much smaller losses than samples sintered fewer times. On the same Fig. 4 the re- 
sults of the measurements of Roberts (1947) made at a frequency of 1 MHz are 
indicated with crosses. The dependence on the temperature has the same chara- 
cter, but the numerical. values are much smaller. Data of other authors lie nearer 
to ours. (Powles 1948, Siniakow 1951, Kazarnovskij 1951). In the case of dielectric 
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losses, there seems that, apart from the technique of preparing the samples, some 
impurities play here a decisive part, as well as the dependence on the current 
frequency and the test field strength. This problem will be the subject of further 
investigations. 

Though the technique of the preparation of samples has a relatively small effect 
on the value of losses, yet it influences very obviously the values of the dielectric 
constant. This is illustrated in Table III, where the results of the measurements of the 
dielectric constant at 20°С (ву) and at the Curie point (e,) for different samples 
are shown. For samples D only the order of magnitude has been given, for its value 
at the Curie point increased so remarkably, that the scale of the apparatus was insuffi- 
cient to make a measurement possible. The values of the dielectric constant for samples 
C have been also included in the Table. 


Table Ш 


As may be seen from the Table above, the values of the dielectric constant for various 
samples differ very considerably, in spite of the fact that the material used for preparing 
the samples was identical. Also during the sintering process no impurities could enter. 
the samples, because the latter, as was already mentioned, were sintered on a liner 
made of the same material. The dielectric constant increases with increasing number 
of sinterings and the time of their duration. At room temperature the dielectric constant 
of samples В is almost 10 times greater than that of samples А and the dielectric con- 
stant of samples D — about 25% greater than that of samples B. At the Curie point 
the dielectric constant of samples B surpasses more than ten times that of samples А, 
and the dielectric constant of samples D almost twice that of samples B. Undoubtedly 
this result is partly due to porosity. At the present moment it is, however, impossible 
to establish whether only porosity plays here a part, or there occur also some other 
changes of the macroscopic structure. Yet the results obtained by us explain in an easy 
manner the discrepancies between the values of the dielectric constants of the ceramic 
BaTiO, met in the works of different authors. | i 

We wish to express here our thanks to all who enabled us to execute the above 
investigations, namety to Professor Trzebiatowski, Professor Roliński., and the 
Kasprzak Institute in Warsaw for providing the necessary materials, to Professor 
Konarzewski and the Fireproof Materials Factory in Skawina for granting us ovens 
for sintering, and to Professor Krupkowski for allowing us to carry out the pressing 
of samples in his Institute. 
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KPATKOE СОДЕРЖАНИЕ 


Ежевский и Пэх, ЗАВИСИМОСТЬ ДИЭЛЕКТРИЧЕСКИХ СВОЙСТВ 
КЕРАМИЧЕСКОГО BaTiO; ОТ ТЕХНОЛОГИИ ИЗГОТОВЛЕНИЯ ОБРАЗЦОВ 
ПРИ ТОКАХ ВЫСОКОЙ ЧАСТОТЫ. 


Исследована зависимость постоянной диэлектрической и потерь керамического 
ВаТ1Оз от способа изготовления образцов. Разработана методика измерения ди- 
электрический постоянный материал с очень большими диэлектрическим 
постоянными при помощи токов высокой частоты, и методика измерения 
потерь. Приготовлено четыре группы образцов BaTiOs, назначенных A, B, С, D. 
Каждая следующая получена была из предыдущей посредством основателнього 
раздробления, сжатия под прессом и вторичной агломерации при температуре 
около 1350? С. Последняя группа D находилась в агломерационной печи особенно 
долго, ибо вместе с нагреванием и охлаждением около 130 часов. Установлено, 
что постоянная диэлектрическая каждой следующей группы образцов, во всех 
температурах, гораздо больше постоянной диэлектрической предыдущей группы. 
Потеры образцов в ряду А, В, С, D, постоянно, но не особенно значительно, 
уменьшаются. 
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ANGULAR DISTRIBUTION OF DEUTERONS 
FROM “Be (p, d) Be. II.“ 


By J. DABROWSKI AND J. SAWICKI 
Institute of Physics, Polish Academy of Science, Warsaw 
(Received May 3, 1955) 


The formula for the Born approximation cross-section for the *Be (p, d) Ве pickup 
reaction obtained with the help of the „Ве + neutron“ model of the 7 Be nucleus is 
compared with the experimental data for 6,5 and 22 MeV protons. The best agreement is 
obtained for the $ Ве — neutron square well radius rọ = 5 · 10713 cm. The angular dis- 
tributions for larger angles lie under the experimental curves mostly due to the com- 
pound state formation. The absolute values of the cross section obtained by taking into 
account the Coulomb correction seem to be too great because of the neglect of the $ Be — 
proton interaction. Comparison with the results of the Butler approximation determin- 
es the limits of validity of this approximation. The same model of the ? Ве nucleus is 
applied for the calculation of the Born-cross section for the ? Be (d, t) *Be reaction. The 
results are compared with experiment for 7,7 MeV deuterons. The “®Be + neutron* 
model of the Be nucleus is discussed. 


Introduction 


The expression for the differential cross section for the pickup reaction ?Be 

(р, d)sBe was obtained in I with the help of the Born approximation and by neglecting 
the 8Be — proton interaction. The complete integrations in the Born approximation 
matrix element leading to the final formula (I, 11) were possible due to the assumption 
of the “Be + neutron“ model of the Be nucleus in P, state with the square well 
8Ве — neutron potential. | | l 
- The aim of the present paper is to compare the results of our calculations with the 


‚ experimental results of Cohen et al. (1953) and Harvey (1951). Thus it is the proof 


of the “8Be + neutron“ model from the point of view of the pickup (p, d) reaction. 
This model, it is well known, satisfactorily explained the electrodisintegration (Mama- 
sachlisov 1943, Guth and Mullin 1949) and photodisintegration at low gamma quanta 
energy (Guth and Mullin 1949). 


* This is the second part of the paper by Dąbrowski and Sawicki (1955a), which will be 
further referred to as IJ. 
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It is well known (Gerjuoy 1953) that the main feature of the Butler (1951) approx- 


imation is to cut off of a sphere of radius го from the г, integration domain іп the Born 
formula!. In order to obtain the Butler expression for the differential cross section the 
first term in eq. (I, 9b) must be omitted the second (multiplied by N.) being retained. 
This method of passing to the Butler expression was more convenient then that 
indicated by Daitch and French (1952, eq. (8)) being also simple in the case of the assum- 
ed square well 3Be-neutron interaction. In the present paper the comparison with 
the results of the Born and Butler theory will be given. For the comparatively great 
proton energy E, — 22 MeV used in experiments by Cohen et al. (1953) it is interesting 
because for the correspondingly large values of the momenta the domain r, Tr, becomes 
important. 

The “8Ве + neutron“ model of the Be nucleus is also aplied to the Be (d, t) Be 
reaction using of the Born approximation. The results obtained are compared with the 
experimental results of Holt and Marsham (1953) and the theoretical curves of Newns 
(1952). 

For both reactions the Coulomb correction is taken into account in the calculation 
of the absolute values of the cross section. 


Numerical Results. Coulomb Corrections 


The ?Be square well potential was assumed to have the radius гу = 5 - 10-13 cm 
(the respective depth being Г, = 12,09 MeV). This value is found most satisfactory 
in the theory of electrodisintegration and photodisintegration (Mamasachlisow 1943, 
Guth and Mullin 1949). Besides the values rọ = 3 · 10-13 cm (И, = 28,92 MeV, 
and го = 7 · 10-13 cm (V, = 7,28 MeV) are discussed. For the deuteron potential 
parameters we take the values / = 21 MeV, а == 2,82 - 10-13 cm. 

In order to compare our results with experiment we obtain the differential cross 
section in the laboratory frame d ode with the help of the known formula (see e.g. 
Schiff 1949): 


do, |+ у? + 2 у cos O, Fi do, 


= 1 
dQ, 1 + у cos O, dQ, @ 
where 
SA ЕСА 2 kp 
у, LAT Va Aid: 9 ka . (la) 


v, is the center of mass velocity and v, is the velocity of the deuteron in the с. m. system. 
The scattering angle O, in the laboratory frame is connected with the c. m. system scat- 
tering angle Os by 


tg OL 26 SW (1b) 


1 The'same symbols will be used as in I. 
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In Fig. 1 the angular distribution of deuterons from the ?Be (р, d) ®Be reaction 
calculated from (I, 11) for ry = 5 - 10-13 cm is compared with the experimental results 
of Cohen et al. (1953) and Harvey (1951). 

In Fig. 2 the comparison of our angular distributions for ry = 3, 5, 7: 10-13 cm 
with the data of Cohen et al. (1953) is given. In Fig. 2 the Butler theory angular distri- 
butions are also given for the same values of 70. 


„% 22 MEV PROTONS (COHEN) 
---- 5-8 MEV PROTONS (HARVEY) 


} CALCULATED 


—22 MEV PROTONS 
65 MEV PROTONS 


(RELATIVE ) 


d6 
10, 


— M 2 
10° 30° 50° 0 E oh 
Fig. 1. Angular distribution of deuterons from the ?Be (p, d) 5Be (ground state) reaction. The experimen- 
tal points are taken from Cohen et al. (1953) 


In Fig. 3 analogous curves are compared to the experimental data of Harvey (1951). 

The eq. (I, 11) makes it possible to calculate the absolute values of the cross 
section. However it is known,that the Coulomb effect lowers these values very definitely 
(Yoccoz 1954, Sawicki 1955). The most important role is played here by the normaliza- 
tion factor of the Coulomb wave functions. The Coulomb effect on the angular distri- 
bution is merely to give some flattening of the curves. In the present paper we shall 
confine ourselves to the estimate of the Coulomb effect on the absolute values of the 


cross section. In eq. (I, II) the plane waves eite p, e were used to represent the 
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incident proton and the scattered deuteron. On taking into account the Coulomb 
effect the plane waves should be replaced by the respective Coulomb functions (see 
e. g. Schiff 1949). The main influence of the Coulomb effect on the absolute values 


of the cross section comes from the factor 


м 27 Np 27 na 

NZ (220) к= (2) 
equal to the product of the ratios of the squares of the Coulomb functions and the 
respective plane waves (at the zero values of the arguments). In eq. (2) n,= 4e? Мы АЗЕ; 


eeee Experimental Points 
=== Bulle 
UY | „ -3 10-90 


( Relative ) 


Fig. 2. Angular distribution of deuterons from the Be (p, d) Ве (ground state) reaction for 22 MeV 
protons. The experimental points are taken from Cohen et al. (1953) 


па = 4e?M;/h®k,. Thus the Coulomb effect will be taken into account by multiplying 
the absolute values of the cross section calculated from (Т, 11) by the Coulomb correc- 
tion factor F. The respective values of this factor are: 
F— 0,132 for E, = 6,5MeV, . 
0,344 for E, — 22 MeV. (3) 
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In Table I the absolute values of the differential cross section for C= 


0 are given 
with F being taken into account. 


Discussion of the Results 


From Fig. | it is seen that the angular distributions for ro = 5 ` 10-13 cm are 
similar to the experimental ones.. This comparison is not quite complete because of 
the lack of experimental data for small scattering angles. However for proton energy 
E, = 22 MeV the theoretical curve falls off too rapidly with increasing scattering 


--- Experiment. 


--- Butler | fo * 310° 


bo 
ay Jo .be 


Born | 55 5 10- 13 cm 


fjem 


Butler 


(Relatıve) 


аб, 
10, 


10° 30° 50° 70° 90º ө, 


Fig. 3. Angular distribution of deuterons from the ®Be(p, d)*Be (ground state) reaction for 6,5 MeV. 
protons. The results of Harvey's experiments for E, = 5, 6, 7, 8 MeV (mean value 6,5 MeV) are taken 
from Cohen et al. (1953). 


angles. It is the common feature of stripping reactions probably resulting from the 
continuous background given by the compound state mechanism. | 

From Fig. 2 И is seen that an increase in го to the value 7 | 02 em gives for 
E, = 22 MeV much too steep a Born approximation angular distribution, and the 
value го = 3 * 10-13 cm gives a somewhat better agreement with the experimental 
angular distribution. 
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On the other hand from Fig. 3 we see, that for E, = 6,5 MeV the value го = 3 
- 10-13 cm gives too flat a Born angular distribution. Conversely the value ry — 
= 7. 10-13 cm yields an agreement similar to that for ry = 5 10-13 cm (however 
for large scattering angles the curve rg = 7 · 10-13 cm falls still faster than does that 
Тогу == 5 ШО em). 

Some flattening of the angular distributions would be obtained by taking into 
account the Coulomb effect (Yoccoz 1954, Sawicki 1955). On the other hand, however, 
the taking into account of the nuclear interaction 8Ве — proton gives some inverse 
effect (Horowitz and Messiah 1953). 

The conclusion is that the value ry = 5 * 10-13 cm gives comparatively the best 
agreement of the theoretical and experimental angular distributions. Thus the simple 
*8Be + neutron“ model of the Be nucleus proposed by Mamasachlisow (1943) and 
then used by Guth and Mullin (1949) receives new confirmation. Apparently the Born 
approximation of the pickup reaction is crude and thus theimportance of this confirma- 
tion is limited. It would be very interesting to compare our calculations with further 
experimental data announced by Reynolds and Standing (1954). 

Let us come now to the comparison of the results of the Born and Butler theories. 
From Fig. 2 it is seen that the Butler theory curves for E, — 22 MeV are in disa- 
greement with the experimental data. This disagreement is the greater the larger is To. 
The fact that the Butler theory is not suitable for the pickup reaction at high energies 
was mentioned by Selove (1955)2, who explains this fact too. Namely the Butler approx- 
imation relies upon the modification of the integral 


fdr, e В (% Dm; (4) 


by the removal of the integration domain r, < rg. We see, that for E, = 22 MeV К 
takes on values from 0,415 - 10!3cm-! for 9, = 0 and for 9, = 60° К = 1,09 - 1013ст-1 
corresponding to the wavelengths 27/k = 15,1- 10-15 cm and 5,6: 10-13 cm respective- 
ly. The wavelengths are of the order of rọ and thus the contribution from г, < rg 
has significant influence on the value of the Fourier transform (4). Thus the removal 
of this domain from (4) gives results disagreeing with experimet. For E, — 6.5 MeV 
the wavelenghts mentioned are equal 24,4 - 10713 cm for O, — 0 and 10,2 - 10-13 cm 


for O, = 60º. Here these wave lengths are greater than ry yet not enough, especially | 


for го = 7 10-13 cm, to justify the removal of the domain г, <r, from integral (4). 
In Fig. 3 it is seen that for E, = 6,5 MeV and r, = 7: 10-13 cm the Butler theory 


curve is in disagreement with the experimental curve. For smaller гу the Butler curves 


become more like the Born curves, but havezeros of angular distribution which do 
not exist in the experimental data (compare Daitch and French 1952). 


Thus the Butler approximation should be treated as a way of avoiding diffuculties | 
resulting from lack of knowledge of the nuclear wave functions in the interior of the 


=—— À——— 


2 The authors are much indebted to Dr W. Selove for letting them see his results prior to р 
cation. 


ee ee 
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— . 


nucleus in the general case. This way is justified only if the internal domain has little 
importance for the values of the Born integrals. It is then not applicable to higher ener- 
gies, for which on the other hand the application of the Born approximation is the more 
justified. 

This fact, visible in the results of our calculations (Fig. 2, Fig. 3), shows the failure 
of the suggestion of Cohen et al. (1953) to compare the experimental data with the 
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0e 8: 
Fig. 4. Angular distribution of tritons from the ?Be(d, 1) Ве (ground state) reaction for 7,7 MeV deuter 


ons. The experimental points of Holt and Marsham and the Butler and Bhatia curves are taken from 
Newns (1952) 


Butler theory. The result given by the Born approximation is the explanation of the 
general resemblance of the angular distributions for E, — 6,5 MeV and for E, — 
— 22 MeV. This result was hard to explain by the Butler theory. 

As concerns the absolute values of the cross sections there is a lack of experimental 
data. One can only guess from the mention by Harvey (1951), that the value of the cross 
section for E, = 5 — 8 MeV is 0,024 barn/steradian for ©, = 20º. The corresponding 
theoretical value for гу = 5 * 10-13 cm is 0,26 barn/steradian i.e. about 11 times too 
great. However it is known from the paper by Horowitz and Messiah (1953), that the 
application of the nuclear 8Ве — proton interaction would lower the absolute value 


i iz 
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of the cross section. In the cases investigated Horowitz and Messiah found a decrease 
by a factor 2 — 6. It is thus hoped that in our case the application of the interaction 
mentioned may lead to agreement between theory and experiment. It is worth while 
to notice, that apparently the Butler formulas give smaller absolute values. These 
values are not discussed here because of the disagreement of the Butler angular distri- 
bution with experiments. 


Table I. Absolute values of (doz /d ) 0, = 0 in barn/steradian for the Be (р, d) Be (ground 
state) reaction. The Coulomb correction is taken into account. Asterisks denote cases in which the 


angular distribution is in distinct disagreement with experiment 


E =6,5 MeV 


3 0,213 0,288 

5 0,134 0,058* 

d 0,054* 0,005* 
Butler Butler 


The %Ве (d, t) Be Reaction 


This reaction was investigated by Newns (1952) with the help of the Butler theory 
and the Born theory with the approximation of Bhatia et al. (1952). Now we investigate 
this reaction with the help of the exact Born approximation with the application of the 


: “Be + neutron“ model of the Be nucleus. 


On neglecting the 8Ве — deuteron interaction and the Coulomb effect we obtain 
for the differential cross section in the c. m. system the final result: 


dos SM, Mal kh. a 
dQ, ^ Qa ы ^ de Te (5) 


where M; = (24/11) M, MZ = (18/11) M are the respective triton and SEMEN reduc- 


ed masses, k, and k, are the respective momenta, 


I = V3 Va fas R (ra) ki, (gra! Vara (6) 
| 6 


where q = (8/9) ki — E, and 


> > > — 2 


3 It is seen from Table I, that the absolute values of the cross section are comparatively little 


sensitive to the values of To, Contrary to the suggestion of our preliminary note (Dabrowski and Sawicki 
1954). 
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where p = k, (2/3) k.; $ = rq — ть и = (V3/2) o. are internal coordinates of the triton; 
X. (u, s) = N,e 27 due) is the internal wave function of the triton used by Newns 
(1952) and y(o) = Мн (ee — e/g is the internal Hulthen wave function of 
the deuteron with the parameter values used e. g. by Malenka et al. (1953). 

The numerical computations were performed for rq = 5- 10-13 cm and 7: 
: 10-13 cm for the incident deuteron energy E, = 7,7 MeV. Unfortunately unlike the 
(р, d) reaction the angular distribution was obtained for гу = 5 - 10-13 cm much too 
flat; for ro = 7 · 10-13 cm the agreement with the experimental data of Holt and Mar- 
sham (1953) is better (see Fig. 4). 


To estimate the absolute values of the cross section similarly as for the (p, d) 
reaction the Coulomb factor 


— | 2лта 27 n, 
pz LE ) p ). (8) 


where n, = 4e? M7/h?k,, n, = 4e? Mey was applied. 


V 


For E, = 7,7 MeV Е = 0,11 and the respective values of = „are 87,12 
mb/ster and 72,30 mb / ster respectively for rọ = 5 ° 10-13 cm and " =7 Uae тте 
The respective total cross sections о, are 110,9 mb and 50,6 mb as compared with 
the experimental value given by Wolfgang and Libby(1952) o“ = 230 mb. Unlike 
the Be (p, d) Be reaction, the absolute values of the cross section obtained are smaller 
than the respective experimental value. In the similar case of the Li (n, t)* He reaction 
the absolute values of the theoretical cross section were also small as compared to the 
common feature of the absolute values of the cross section in stripping (Dabrowski, 
Sawicki 1955). 


Remarks Concerning the Model of the 'Be Nucleus 


It is well known that in general the shell (two body) model of the *Be nucleus is 
very crude and does not explain many processes and properties of this nucleus (e. g. 
its magnetic moment). For a description of them a more accurate model “2a + п“ 
ought to be used. 

However the mechanism of the Be (p, d) Be and Be (d, t) ®Be reactions like 
the mechanism of the Be (y, n) ®Be reaction pretty well distinguishes the external 
neutron and indicates the use of the shell model. The fact that the external neutron 
in the Be nucleus is loosely bound with the energy 1,63 MeV compared with 8 — 10 
MeV average binding energy per particle, is visualised in the large values of the 
cross sections for the (p, d) and (d, t) reactions. 

However it is necessary to mention, that recently the “Be-+neutron“ model was 
questioned in the case of photodisintegration even for y — quanta energy 6 MeV by 
Carver et al. (1954). Their experiments suggest, that only 20% of the Be nuclei from 
the 9Be (y, n) Be reaction are produced in the ground state. In the remaining cases 
appear ®Be (in an excited state) and He in the ratio 1,2:1. Apparently in both the 
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(р, d) and (d, t) reactions the competing mechanism does not interfere with the pro- 
duction of ?He. Similarly the reaction leading to ®Be* is eliminated by the proper selec- 
tion of scattered particle energy (e. g. Cohen 1953, Holt and Marsham 1953). 

Further the theory of the (y, n) reaction with the use of the “Ве + neutron" 
model given by Überall (1953) for large у — quanta energy 20 MeV < E, < 200 MeV 
was found to be in complete disagreement with experiment (Jones and Terwilliger 
1953). Similarly the “8Ве + neutron* model may be less important for large energies 
E, and Е. However for the energies E, = 6,5; 22 MeV and EA = 7,7 MeV the agree- 
ment with experiment was found satisfactory for rj = 5 — 7: 107? cm. The best 
value seems to be rọ = о. 10-13 cm.*. 

We are much indebted to Prof. L. Infeld and Prof. W. Rubinowicz for a helpful 
discussion and their kind interest in this work, and to Mr. W. Czyz for interesting 
discussions. 


APPENDIX 


On the Born Formula for the Cross Section 


We shall give an outline of the deduction of the Born cross section in the с. m. 
system for the ?Be (p, d) ®Be reaction in order to explain the simplifying assumptions 
used in this deduction. 

The hamiltonian of the system Н may be written down in two forms 


Н= H,(£r,o,) + Г, + Tz, = Н.) + Ha ( c, o) + + Ti (1А) 


Ho, Hs, H, are respectively the internal Hamiltonians of the Be, Ве nucleus and 
the deuteron. T7,, Tr, are respectively the kinetic energy operators of the systems: 
proton — Be nucleus and deuteron — 9Be nucleus. Vz, V, are the proton — neutron 
and Ве — neutron interaction potentials. In (LA) the Be — interaction potential is 
neglected. 

We try to solve the Schrödinger equation of the problem 


H y — Ey (2A) 
in the form | 
V = ж T Vi. (3A) 


where yy is the plane wave of the proton incident on Be and y, represents the deuteron 
departing from the 8Be nucleus. 
On substituting (3A) in (2A) we find with respect to (1A): 


(Hs + Hi + T + Vo — E) v1 = — Vay. (4A) 


— — —ñ æ —E＋— q —ͤ —„- ß 


* Although the value го = 7: 10-13 cm gives a better agreement with experiment than the value 
ro = 5 * 10718 cm for the angular distribution for the (d, t) reaction, the value gives a better value for 
the total cross section for the same reaction. 
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The potential V4 is contained in H, and thus on the left hand side of (4A) we have 
(Va + Vo) yı. The first Born approximation rejects, the term? уру. Thus we neglect 
the term Voy, as compared to V, y,. This corresponds to the mechanism of the (р, d) 
reaction, in which the neutron is picked up from the Be nucleus due to the У, potential 
and in spite of the Йу potential. Besides the n — p interaction is stronger than the n — 
Be (greater binding energy). It seems that taking the term Йу y, into account would 
lower the cross section. It is however beyond the first Born approximation. 

On neglecting Г, in (4A) we can calculate y, in the usual way and then the cross 
section. 

It is interesting, that in the calculation of the cross section for the inverse 8Ве 
(d, p) Be reaction V, is neglected as compared with V, in the equation analogous to 
(4A). Notwithstanding this both the cross sections are connected with the “reciprocity 
theorem for nuclear reactions*. The validity of this theorem for the Born approxi- 
mation can easily be proved explicitly (Bransden 1954). 


KPATKOE СОДЕРЖАНИЕ 


Домбровский и Савицкий. УГЛОВОЕ РАСПРЕДЕЛЕНИЕ ДЕЙТРО- 
НОВ С РЕАКЦИИ 9Ве(р,а)8Ве II. 


Полученный при помощи модели ,‚8Ве+фнейтрон” формулы эффективного 
сечения в приближении Борна для реакции „pickup“ ?Be(p,d)PBe сравнено 
с опытными данными для протонов с энергией 6,5 и 22 MeV. Самую лучшую 
согласованность получено для прямоугольной ямы потенциала „#Ве+нейтрон” 
с радиусом ro=5.10-13 см. Угловое распределение для больших углов лежит под 
экспериментальными кривыми главным образом благодаря образованию сложе- 
ных состояний. Абсолютные значения эффективного сечения, полученные учи-. 
тывая куломбовскую коррекцию, кажутся чересчур большими по причине упу- 
щения из виду взаимодействия Be и протона. Сравнение с результатами, полу- 
ченными на основании приближения Батлера, определяет предел применимости 
этого приближения. 

Эту же самую модель ядра Be применено при вычислении борновского 
эффективного сечения реакции ?Be(d,t)!Be. Результаты сравнено с опытом для 
дейтронов с энергией 7,7 MeV. Придискутирована модель ядра Be ,,#Ве--нейтрон”. 
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ELECTRONIC INTERACTION IN THE FREE-ELECTRON MODEL 
FOR NON-BRANCHED CONJUGATED SYSTEMS 


Ву S. OLSZEWSKI 
Institute of Physics, Polish Academy of Sciences, Warsaw 
(Received May 30, 1955) 


The ASMO CI method for molecular orbitals of a free-electron for non-branched 
bond systems of conjugated molecules is presented here in close analogy to the ASMO CI 
theory as developed by Goeppert-Mayer and Sklar and Craig for molecular orbitals in the 
form of the linear combination of atomic 2p orbitals. This scheme was used to calculate 
the energy levels in two cases: ethylene and butadiene. The results are in a better agree- 
ment with experimental data than those which can be obtained by analogical computations 
with LCAO MO method. This is the case especially in the singlet-triplet separation for the 
first excited states. In addition a significant decrease in the configurational interaction 
effect compared with the LCAO MO method was obtained. 4 


1. The Theorv 


The theorv of electronic spectra of unsaturated hvdrocarbons in molecular 
orbitals (MO) formulation-is generallv used in two wavs: the first is the linear combi- 
nation of atomic orbitals (LCAO), (see e.g. Pullman and Pullman 1952), where the 
molecular orbitals are built from the atomic 2p orbitals of individual atoms. In the 
second one, a free electron (FE) model, we have the quantization problem of free elec- 
tronic motion in a long thin tube extended along the molecular free-bond path. In 
these two methods the - electrons are treated independently from the rest of the mole- 
cule which manifests itself only as an effective “core“. It is in the field of this core 
that the z-electrons move. (Ruedenberg and Scherr 1953). 

In the LCAO MO method in its simplest form the presence of the core is taken into 
account in the Coulomb and resonance integrals usuallv extrapolated from the exper- 
imental data. The FEMO method in the form as it is at present developed assumes the 
interaction potential of the core with z-electrons as constant. This potential may be 
taken as equal to zero. The only quantity taken from experiment in the latter method 
(and in the exact LCAO MO method as well) is the mean length of the bond between 
the carbon atoms in the conjugated systems of double bond.! 


1 Recently there appeared some papers which confirm a close correspondence between both FE and 
LCAO treatments (see e. g. Ruedenberg 1954). These results may be considered as a great support for 
our theory. 

(419) 


The next step in the LCAO MO method isto build the wave functions correspond- 
ing to the particular electronic configurations of a molecule. These wave-functions 
take the form of antysymmetrized products of molecular spin-orbitals (ASMO). 
At the same time in the Hamiltonian operator mutual electronic interaction terms 
are included. (Goeppert Mayer and Sklar 1938). 

In order to determine the energy of the stationary states of the molecule as 
precisely as possible the configurational interaction (CI) is taken into account. The 
full л-е1есїтопїс wave function of a molecule is built for a given state in the form of 
a linear combination of the eigenfunctions of the configurations, (Craig 1950). 

In the calcultations all integrals may be computed by assuming only the С— С 
distance as known. Thus a purely theoretical scheme of a consistent theory of electron- 
ic spectra for the molecules investigated is obtained. We shall call it further on the 
“exact LCAO MO method“. 

The aim of this work is to develop the theory of electronic spectra of unsaturated 
hydrocarbons by means of a free-electron model in a similar way to that of the con- 
ventional LCAO MO method. We shall carry it out first for a non-branched model. 
The generalization of the theory for a model describing branched bond path systems 
is in preparation. 

The purpose of the free-electron theory is to solve the eigenproblem in the 
following form 

Hp = Ey (1) 
2 


h 
where H= ә, V^ --V- 


which corresponds to the electronic motion in a box extended along the molecular 
core in rectangular coordinates. 

We assume, as usual, that the potential is only dependent on the x bond path 
coordinate corresponding to the core, i.e. 


V = Vix) | 


We assume further the size of the cross-section of the box to be small in comparison 
with its length. This means that the excitations energies in the y and z directions are 
much higher than in the x direction. Thus the movements in y and z remain in the 
ground state while that in x direction becomes excited. (Ruedenberg and Scherr 1953). 


By this assumptions our problem has the three-dimensional solution: (Rueden- Ч 
berg and Scherr 1953). 


p (% у, 2) = p(x) (Je) sin (лу/е) sin (2mz/e), (2) 


where г is the size of the square cross-section of the box, and the following relations 


are fulfilled 
0 X y Se, — le Xz e, e<l 
x as before being the coordinate tangent to the path bond, y being normal to x in the 
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plane of the molecule, and z being normal to the plane of the molecule, which we take 
as z — 0. v represents the one-dimensional molecular orbital depending on x. 
For a molecule containing s z-electrons moving along the core, the Hamilton oper- 


ator takes the form 
1 
Qe == 3 Hu Es 2 Y ear (3) 
u 


wy 
where H“ is the Hamiltonian (1) for the u-th electron moving in the field of the core 
itself. 

For the non-branched molecular model the z-electron motion along the broken 
line of the core (e.g. in polyens) can be approximated by the motion in a straight- 
linear box with a length equal to the free bond path length (see Ruedenberg and 
Scherr 1953, Fig. 2). Then we have 


Ги» = ү (x, bn 3:5 + (Yu = у)? + (2, E 2,)° (4) 
The perturbations introduced by this assumption may be taken as insignificant and тау 
be neglected (particularly for the all-trans case). 
The wave functions of the s electron system are built from normalized antisymme- 
trized product functions of the type (Mulliken 1949) 


(фу а)! (pı B)! (p; o)! .. | 
Da = (6% | (ру a)? (ру B)? (a c).. |, (5) 


sesssssss cct] n 


where 4 is the running index by which the correspondence of the electrons to the 
particular molecular orbital q, together with the spin function a or f is characterized. 
In other words, 4 denotes the spin-orbital configuration. 

The energy of the electronic state corresponding to the wave function D, takes 


the form (Mulliken 1949) | 
1 ; 
ға= Yo +, Vim — Kw) (6) 


(the summations runs over all occupied spin-orbitals), where 


I, = f p, (u) Н“ o, (и) dr (2) 


is the core energy for molecular orbital , equal to the corresponding eigenenergy 
for the unperturbed problem, 


Jam = | 0: ж 0) Eos Pm 0) a, Q 


being the Coulomb integral between the molecular orbitals 9, and q, and 


Kim 2, ji Pn (и) E! (v) L Pm (и) Pn (v) dt, йт, (9) 
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being the exchange integral between g, and Pp, where K, = Kum when q, and 
q,, have the same spin function in Ф and K, = 0 when they have different spin 
m 5 


functions in D4. as 
Now, as the length of the cloud of the 7 electrons is large compared to its width, 


we may assume that approximately 


— x. (4a) 


In view of the above considerations and the normalization of the ꝙ term depen- 
dent on y and z we get from (7), (8), and (9) 


I, = f v; G HY v, (e) dx, — Em (та) 


where E, is the eigenvalue corresponding to the Ф, orbital, 


sj -fv (x) Va (x,) | 


т» = |X, 


2 
1 Vn (Xu) Vm (х,) dx, dx, (8a) 


Хи 
апа 
* * e? 
Кы = ji Wa (Xu) Ym (№) Кү" (хи) Yn (v) day dx. (9a) 
Xu — № | 
In the case of a non-branched bond path we have the following boundary conditions 


y (x) = 0 (10) 
on both ends of the path. 
From (10) and from the normalizing condition 


Д 
[$e d 1 (11) 
0 


we get the following set of eigenvalues and eigenfunctions of the unperturbed problem 
(V(x) = 0) 
h? n? 2 


E Setran ee, 0408) 


We assume that the free electron path exceeds for a distance C — C beyond the last 
carbon atom in the conjugated system. The formulation (10) — (13) was first given 
by Kuhn (1949). 

Now we shall evaluate the Coulomb and exchange integrals. We shall begin with 
the K integral. For molecular orbitals we assume generally normalized and orthogo- 
nalized functions of the type 


(REI 


Vy; = А; sin ах where O Xxx (14) 


We must first find the most convenient form for the expression IPA — x. We 
get namely? | 


2 This transformation may be easily deduced by means of 6 function theory. See e. g. Ivanienko 
and Sokolov, 1953. 


Electronic Interaction in the Free-Electron Model 495 


со 2л 
1 = 1 | ek | Xu — Xy | 
Er = Gp 2 д? 3 k2 -+ r2 dk rdr dO (15) 
—co 0 0 


Now we can evaluate our K; exchange integral by putting for the molecular orbitals 


of a free electron (14) into the expression for Ki. Thus we have 
ist 


2 
= : : € $ . 
pes "n sin a; x, A;sin a; x, 1 A; sin a; x, А; sin a; x, dx, dx, 
B 0% H 5402 
where x, and x, denote the coordinates of the first and second electron respectively. 
Integration is performed along the whole free electron path. Denoting a; + aj by kı, 
and а; — a; by №, we get for the last formula the form 


11 
L AA 
ETT A; А; [cos kı x, cos ky x, — cos Б x, cos ky x, — 
00 


72 
-— dx, dx, (16) 


| x Ton = № | 


Each component of this integral may be calculated separately. We shall consider 


— cos k x, cos ky x, + cos ka x, COS ka 2,] 


first the integral 
11 


e 
= 3 . ' 17 
C [feos kı x, cos Ё, = puerum dx qas (17) 


00 
which by formula (15) o 


+0000 


eik Ix, x 
C= — cos ax, cos cx, | | oux TA dkrdrdx, dx, 


For clarity we have E k, by a, and T by c. After performing the succesive 
integrations, we obtain the following expression | 


л л? c+ а с si (le) cos lc — а si (la) cos la 
Е C= IE = cos? 3 I + л [cos cl + cos an) il cos ا‎ a la) cos (18) 


In the course of integration we have taken into account the boundary conditions 
for the free electron model (10). It may easily be seen that the result is symmetrical 
in relation to c and a. 

To calculate the exchange integral K, we must also find the value of (ze?) C 
when c--a. This value will be equal to the limit of the last expression when c — a. 


оо+оо 11 


ПШ е ＋ eit 1x, x cos ax, cos ax, rdr dk dx, dx, 
co— oo 0 
л? 1 ШЕН 
cem cos? al + 2л - al [t si (lt) cos lt] t= a- 9 1 al · Si (la) (19) 
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Evaluating the derivative, we have also considered the boundary conditions for our 
problem. We get obviously a similar expression when we have c instead of a. 
Taking advantage of the above results the Coulomb integral J for our problem 
will easily be evaluated. In our case it has the form 
Ll 
e? i 
Ji = ffa sin a; x, А; sin aj 2, ° T EU A; sin a; x, A; sin а; x, dx, dx, 


00 
1 


1 
= Zh 4 ffa — cos 2a; x, — cos 2a, x, + cos 2a; хи * соз 2a; ,) dx, dx, (20) 
00 


The last integral is quite similar to that evaluated in (18); the second and third ones 
readily follow from it. The first element in Jj, i. e. the integral 
оо+оо 1 1 


NE dr dk dx, d 
— T ar X. X, 
po 2 т 
—оо 00 nm 


0 


may be evaluated by finding the limit of the expression 7 cos? al Si (la) when . 
a 


The result is xl. The evaluation of the % type integral presents no difficulty. 
The final form of the integrals in the 45 FEMO method for the non-branched 


model is 

ef cos? ki «. cos? ko | <. 2x 2 it he 
n E Si (i) + US Si (lko) EE: cos 5 
k; si (lk,) cos lk, — Е, Si (I) cos 1 + 


ki — № 


l 


Kj 4} А}. 


— 2 [cos k, l + cos k,l] ° (21) 


e? л 1 
7) | + Di cos? k,l + 2 cos kal cos kal x 
: (ky + ka) si (kı + А) Lcos (kı + ka) L— (kı — k) si (kı — ko) l cos (ky — ko) l 
ky q ka 


Ji = A; A; 2 


cos? (ki + K E 
_ hh 


cos? (k, — kg) 2 
[a 2 cos (k — №) lsi (k — А) Д (22) 


[ж + 2 cos (kı + ka) 1 si (E, + ka) 1] 


ае? cos? (20,1) c. cos? a;l 
Ja = A: E | + “ae Si (2a;l) — E [z:-I-2cos 2a; l si (2a; 2) (23) 
k, and k, have the same meaning as previously described; а; is the coefficient at x 
for a given molecular orbital. 


Finally we must consider the electronic configurational interaction problem in AS 
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FEMO method. The singular ©, wave function may be, or may not be, а good approx- 
imation to the complete z-electronic wave function О; of the given state of the 
electron system in the molecule. When this function suffices, the energy may be 
computed directly from Eg. (6), but when several О; have to be taken into account, 
we must include the configurational interaction by applying the variational calculus 
to a linear combinations of Ф, functions. 

The О; eigenfunction of the state may be written in the form 


9 = У 4%, (24) 


where Ф; denote AS FEMO wave functions for a given type of symmetry, e. g. suc- 
cessively with the increase in energy. The energies E, correspond to G.. Now if we 
introduce the interaction between the configurational states and if we have as many AS 
FEMO wave functions as possible, then the best energies and wave functions will 
result, this being consistent with the primary application of the free electron molec- 
ular orbitals. 

Thus, in order to obtain the eigenvalues for the energy E and the corresponding 
eigenfunctions, we must solve the secular equation 


EI - E (CT Or) (Ф, | u 


© | 2) Ea a cce ONES (25) 
where (Dr | Фи). = [| I fi y: an dr, dr,. 
u<v 


The entire method is strictly parallel to that given by Craig (1950) for antisymme- 
trized molecular orbitals LCAO. 

Now we ought to evaluate the non-diagonal matrix elements of our secular equa- 
tion. They are built of integrals over molecular orbitals of the type 


[j| kl] = [| y; (u) gj ш= qx (0) pı 0) dz, dr, (26) 


For the non-brached free electron model, we get for this integral the particular form 
2 
с е . А 
A; Aj Ar Ai lil sin а; x, sin a; x, EREN sin Ak x, sin aj x, dx, dx,. 


From the above it may readily be seen that this integral may be expressed in terms 
of the integrals computed previously. These integrals will be of the C type. 


2. Two simple examples: ethylene and butadiene 


Though the AS FEMO method may be easily applied to much larger molecules, 
the ethylene and butadiene molecules have been chosen as examples not so much 
with regard to their simplicity as to the highly developed experimental investigation 
of their spectra and precise calculations by the LCAO method. 
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The entire method developed above has been applied to both ethylene and 
butadiene. The results have been compared then with those obtained by the far more 
complicated non-empirical method of antisymmetrized combination Slater 2pz 
atomic orbital products and experimental data. In the computations only the free 
electron molecular orbitals corresponding to the LCAO MO have been considered. 

For ethylene we take two free electron molecular orbitals of the form 


V Ver т 47 
| = jiet Ex 
5 en 
y r 
Vo = 55 


Hence we get the energies of the configurational states (Parr and Crawford 1948). 


1 = 2h + Ju 

2 = 21,4 Jo 2 
28 

Vie 9 


АЛ Ki 


12 


(1 | 2) = K, is the nondiagonal element in the configurational interaction matrix, 
I, and J, are taken from (12), and Jj, Joo, Лә, Кә from (21) — (23). 

The numerical values are given in tables I —III. The 1 and 2 states have the con- 
ventional meaning: we assign two electrons to y, and wy, orbitals respectively. Т is 
the triplet and Г the singlet state which are obtained when there is one electron on 
each y, and ys. N and Z states are obtained when the interaction between 1 and 2 
states is taken into account. 


Table I. The electronic interaction integrals for ethylene (in eV) 


Ja = 1,14779 12 = 1,69665 
Jon = 2,28848 Ki = 0,93242 


Table II. The electronic energy levels for ethylene (in eV) (without configurational interaction) 


exact LCAO MO method? 


(Parr and Crawford 1948) FEMO method 


0 0 
12,5 14,90 
1,8 6,40 
10,2 8,26 
—1. 3 — 0,058 


The interaction of 1 and 2 states in ethylene gives a decrease of the ground state 
only by about 0,06 eV, while the configurational interaction of the exact LCAO 
method gives for a similar decresae 1,3 eV. 


SEE EURER DA TE Чү! e DO TE жа Мокін Та) 
$ The values are calculated from Slater atomic orbitals with the effective charge — 3,18. 
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Table III. The electronic energy levels for ethylene (in eV) (configurational interaction included) 


exact LCAO -MO method! 
FEMO : | 1 
(Parr & Crawford 1948) оң We UE. 


0 0 0 

3,1 6,46 3,1—5,64 
11,5 8,32 7,6 
15,0 15,02 


The length of the free electron path in ethylene is taken to be equal to the 3-fold 
double bond length of this molecule, i. e. / = 4,05 A. | 
For butadiene we take 4 free electron molecular orbitals (with successively in- 


creasing energies) 
2 
Vn = Vie 535 (29) 


The lowest singlet configurational state at the energy Eo with wave function V, 
has two electrons assigned to both y, and ws. Excited configurational states arise 
when the electron is raised from y; to h. We shall denote such states by V;and T; 
corresponding to singlet and triplet states respectively. We shall restrict ourselves 
to one-electron excited states. On the basis of (6) we get the following expressions for 
the energy differences 


Е — EO = 15 — I — Jog — 2 Л» + 2h + Љ + Ki — Kis + Ka; 


23 


z frl — EO = I- H Л: > 2e + 2/4 + Ла + Къы — Ku Кү, 
14 (30) 
a — Ey = I, — В — J — 2h + 2Ла + Ла + Кы — Ku £ Ка 
24 
ا‎ — Eo = 1—1, — Л, — Л» + 25% + Jia + Kia — Кз SE Куз 
13 
where the plus and minus signs correspond to singlet and triplet states respectively. 
The energies obtained for the states corresponding to different assignment of 
electrons to particular molecular orbitals may now be corrected by taking into account 
configurational interaction. As we are interested here in the lower excited states, 
only the lowest configurations of a given symmetry have been considered. It seems 
that by including a large number of configurations we would get only a small cor- 
rection of our results (see the discussion below). 
So we take for „mixing“ Vo, V4, and Jia wave functions, all three of 14 symmetry, 


4 The location of the lowest triplet state in ethylene is not quite clear yet; its symmetry is also not 
very certain. The strong absorption in ethylene at 7,68 eV and the symmetry of this state is undeniable. 
The investigations in the liquíd nitrogen temperature (see Moser 1953) show the appearance of the triplet 
state at a relatively small distance from the singlet state — about 1,00 — 1,25 eV (Potts). 
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V and Vj, of 1B symmetry, Тз and Ту of ЗВ symmetry, and Тиз and 124 of ЗА 


symmetry”. 
Table IV. The electronic interaction integrals for butadiene (eV) 

Ju = 0,66408 Jiz = 0,98163 Ja = 1,44256 Jog = 1,48901 

Jos = 1,32405 Ju = 1,08311 hs = 1,05366 

Ky, = 0,53947 K» = 0,73701 Ky, = 0,064591 

Kis = 0,024884 Ки = 0,11659 
[11 | 13] = 0,31755 [11 | 24] = 0,38957 [13 | 24] = 0,07136 
[22 | 24] = 0,46093 [12 | 23] = 0,60594 [23 | 14] = 0,13106 
[22 | 13] = 0,34244 [12 | 14] = 0,06647 [12 | 34] = 0,64098 


The approximate values for the matrix elements of configurational interaction (in eV) 


(Vo | Vis) = 0,57 (Vis | Va) = — 0,50 
(Vo | Voa) EE 1,66 (Vos | Via) = — 0,38 


Table V. The electronic energy levels for butadiene (eV) (without configurational interaction) 


Confi ‘onal exact LCAO MO method б 
onfigurational state (Coa od irese 151) FEMO method 


5 The same system of electronic configurations in butadiene was considered by Pariser & Parr 
(1953) in their semi-empirical method. The experimental results for both molecules are taken from 


this paper. 


(Tis | T;,) = — 0,64 
(Tos | Tia) = — 0,64 
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The numerical results for configurational interaction matrix elements as well as 
electronic repulsion integrals and energies of the states are given in the tables below. 
By solving four secular equations corresponding to the four types of symmetry, we 
get the final results for energies of the electronic states which are given in table VI. 

The length of a free electron path in butadiene was taken to be equal to the 


five-fold mean length of the bond in the conjugated molecules, i. e. 7,0 À. Substitut- 
ing the parameters of the molecular (14) — type orbitals, it may be shown that for 
these orbitals all the values connected with the electronic interaction in the corre- 
sponding molecules (exchange integrals, Coulomb integrals and configurational inter- 
action matrix elements) are inversely proportional to the length of the free electron 
path$. 

By the method given above any molecules having nonbranched free electron 
bond paths may be calculated’. No other assumption concerning the empirical quanti- 
ties than the length of the bond need to be introduced. A generalization of the theory 
for any unsaturated hydrocarbons with conjugated systems of double bonds will be 
presented later. 


3. Discussion 


The first and the most important drawback of the only purely theoretical 
method, developed so far for investigating the spectra of unsaturated conjugated 
hydrocarbons, is the disagreement of the results with experiment. Few molecules 
have been examined so far by this method and the results have been compared with 
the experimental data which are far from being complete. Full computations have 
been made only for ethylene, butadiene, benzene. In spite of their incompleteness 
the present experimental data point to serious drawbacks in the LCAO method. 

In the present work which was concerned with the molecules containing un- 

. branched chains of conjugated bond systems we discuss only the structure of the 
absorption spectra for ethylene and butadiene, where a more exact comparison of 
the given data with LCAO method and experiment is possible. 

The energy levels for ethylene theoretically found by the LCAO method are 
qualitatively in agreement with experiment but quantitatively they are far from that. 
By AS FEMO method the error in estimation of the first singlet excited level for 
ethylene is reduced from 3,9 eV to about 0,7 eV. 


6 The free electron model of butadiene corresponds rather to the trans form of the molecule. The 
cis form could be considered as characterized by a somewhat shorter free electron path. Then the energy 
of the ground state and the energies of the excited states of the trans form would be lower than the corres- 
ponding cis energies, since the “core“ part in the expressions for the energies and the part of the electronic 
interaction are positive and increase with decreasing l. 

7 The free-electron treatment of polyens (with V — 0 assumption) is not quite permissible. 
But for the case of the short’ chain the spectrum is not greatly effected by the localization of the 
a-electrons on the double bonds and this effect may be neglected. 
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It is well known that for butadiene, although the interaction between a large 
number of configurations is included in the calculations with LCAO method (Coulson 
and Jacobs 1951), even a qualitatively correct scheme cannot be obtained. The low- 
est excited singlet state is there a 1A state obtained from the „perturbed“ electron 
jump from level 2 to 4. Vis and Voa, when the configurational interaction is included, 
give states, the difference between them being as high as 7 eV. The results of the 
calculations by the self-consistent method, (Parr and Mulliken 1950) taking into 
account to a great extent the effect of configurational interaction, give at the same 
time the reverse succesion of these states. Furthermore, the lowest singlet excited 
state is in butadiene that of 1B symmetry lying onthe LCAO scheme just over 3,5 eV 
above the first excited 1A state. 

With such an unintelligible description of the positions of the energy states, 
a good quantitative agreement cannot be expected. The first excited state observed 
lies 1,2 eV lower than the lowest 1B level of the LCAO method. For the second state 
in this method a corresponding term can hardly be found. 

Now the present method gives for butadiene a reasonable succession of states in 
agreement with the results of the semi-empirical Pariser and Parr theory (1953) 
and with the calculations performed by means of a self-consistent field method by 
Parr and Mulliken (1950). 

As far as the agreement with experiment is concerned, it seems to be even better 
than that for ethylene. The difference between the results of observations and those 
obtained from computations is about 0,4 eV, with only a small number of configura- 
tions included in the calculations. 

The separation between the singlet and triplet level present another problem. In 
the LCAO method they are far too large. This shows this method to be unreliable 
in evaluating the exchange integrals. The real singlet-triplet separation for the lowest 
excited states is probably approximately eqnal to 1eV (Parr and Mulliken 1950). 
This result is in no great discrepancy with the present method. 

The characteristic feature of the method is that it gives relatively small values 
for the electronic repulsion molecular integrals. Thus the corresponding matrix ele- 
ments, when the configurational interaction is taken into account, have smaller val- 
ues than those built from molecular orbitals, which are the combinations of atomic 
2p orbitals. This results in a descrease in the electron configuration interaction effect 
on the value of the energies of the states. At the same time the configurational inter- 
action method has here a more perturbational character than for LCAO orbitals. 
The decrease of the ground state caused by the configurational interaction for the cal- 
culated molecules is many times less than that evaluated by the LCAO method. This 
fact allows at once for some insight into the system of energy states of the spec- 
trum. The probability of mixing the states given by the theory is strongly reduced. 
Thus the necessity of applying a complicated mathematical formalism connected 
with the interaction of many configurations is also decreased. 

From the above considerations it may be supposed that in our examples the 
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LCAO method based on the calculations of 2p atomic orbitals is not so well confirmed 


by experiment as the FEMO method delocalizing completely the л electrons along 
the free bond path of the molecules. 


The second serious obstacle in applying the purely theoretical LCAO method 
lies in its computational difficulties. The present method avoids this difficulty to a 
great extent. АП the integrals are computed in the molecular orbitals. Thus the neces- 
sity of expanding into integrals over atomic orbitals is avoided. There are no multi- 


centered integrals, and for the integrals applied in the present method simple general 
formulae can be obtained. 


The author wishes to express his deep gratitude to Dr W. Kolos for assistance 
in preparing the manuscript and many stimulating discussions. 


KPATKOE COJIEPZKAHME 


Ольшевский. ЭЛЕКТРОННОЕ ВОЗДЕЙСТИЕ B МОДЕЛИ СВОБОДНОГО 
ЭЛЕКТРОНА ДЛЯ НЕРАЗВЕТВЛЕННЫХ СОПРЯЖЕННЫХ СИСТЕМ. 


Представлен метод АЗМО СТ для молекулярных орбит свободного электрона 
в случае молекул о неразветвлённой системе сопряженных связи, при полной 
аналогии с теорией АЗМО СТ, развитой Гоепперт-Майер и Скляром и Крейгом 
для молекулярных орбит, как комбинацией линейных атомных орбит 2р. 

Эта схема была применена к вычислению энергетических уровней в двух 
случаях: этилена и бутадиена. Полученные результаты сходятся с эксперимен- 
тальными данными более, нежели результаты, полученные при помощи анало- 
гичных вычислений методом [САО MO. Происходит это в особенности при сепа- 
рации синглетно-триплетной первых возбуждённых состояний. Кроме того полу- 
чено значительное уменьшение эффекта конфигурационного возддействия в срав- 
нении с методом LCAO MO. 
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ЕВКАТА 
Acta Physica Polonica, 14, 407 (1955): 


ANGULAR DISTRIBUTION OF DEUTERONS FROM “Be (p,d) Be 


Bv J. DABROWSKI AND J. SAWICKI 
Institute of Physics, Polish Academy of Science, Warsaw 


The angular distributions of deuterons from the ®Be(p,d) Ве reactions and tritons 
from the ?Be(d,t)$Be reaction were analyzed in our previous papers (1954, 1955a, 
1955b, 1956). The results of these papers require a correction. Namely, the expression 
for L(k) adopted from the paper by H. Überall contains an error overlooked in our 
previous papers. In Eq. (I, 3b), in the denominator of the expression for N?, n, should 
be replaced by пф, and the opposite sign should be adopted for Ni. The sign of the 
term k sin kr, (1 -- Нато) in Eq. (I. 9b) should be changed to positive. 

The most important correction concerns the absolute values of the cross sections. 


Table I of paper (II) and Table II of the Supplement (1956) should read: 


TABLE I 


E, — 22 MeV 


Butler Butler 


TABLE II 


Theoretical values 
for ro = 5 X 1013 cm 


For E, = 5—8 Me V and 0, = 20°forry = 5 X 10-13 cm (Born) one ARE do,/d2, = 
= 0. 22 barn/ster, i. e., a vate about 9 times too great — instead of 0.26 barn/ster, 
i. e., a value about Il times too great. 

For the 'Be(d,t) Be reaction (doç/dQO s)oç=00 = 34 mb/ster and 39 mb / ster 
instead of 87.12 and 72.30 mb / ster for гу = 5 x 10-13 cm and rp = 7 x 107? cm respec- 
tively. The respective total cross sections are 44 mb and 27 mb (instead of 110.9 mb 
and 50.6 bm). The angular distribution for this reaction remains unchanged. 
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The agreement of the calculated angular distributions with the experimental 
data for the (p, d) reaction is only very slightly changed (see the corrected Figs 2, 3, 5). 
The only difference is that the best value of ro should be somewhat smaller than 
9 X 1013 cm. The curves for г = 7 x 10-13 cm are omittedi n the corrected figures 
since they are in distinet daisgreement with experiment. 
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Fig. 5 (corrected) 
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